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X> : ABSTRACT 

Centaurus X-3 is a well-studied high-mass accreting X-ray binary and a variable 
source of high energy gamma rays with energies from 100 MeV to 1 TeV. Previous 
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results have suggested that the origin of the gamma rays is not the immediate vicinity 
of the neutron star but is sited in the accretion disc, perhaps in an accretion wake. The 



Durham Mark 6 gamma ray telescope has been used to measure the gamma ray flux from 
Centaurus X-3 with much higher sensitivity than previous ground-based measurements. 
The flux above ~ 400 GeV was measured to be 2 ± 0.3 x 10~ n cm~ 2 s _1 and appears 
constant over a period of 2 - 3 months. In 10 hours of observations there is no evidence 
for periodicity in the detected gamma rays at the X-ray spin period either from a site 



Q H i in the region of the neutron star, or from any other potential site in the orbit. 
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The high-mass X-ray binary Centaurus X-3 was the first binary from which coherent X-ray 
pulses were detected flGiacconi et al. 1971 ). It contains a 4.8 s pulsar in a 2.1 day orbit around 



an O-type supergiant, V779 Cen. The pulsar period has been shortening since discovery; this is 
attributed to accretion of matter by the neutron star from the more rapidly rotating inner edge 
of its accretion disc. The long-term period history was monitored by GINGA and shown to be a 



linear variation with a small added sinusoidal term of about 9 year period ( [Tsunemi et al. 1996[) 



However, this long-term periodicity has not been observed in recent data from BATSE ( Bildsten et 



al. 1997 ). Large short-term variations in both X-ray amplitude and instantaneous period derivative 



are seen with no apparent correlations ( [Finger, Wilson fc Fishman 1995 ; Tsunemi et al. 1996) 



suggesting that the X-ray luminosity observed is often only a small fraction of the accretion power. 

Very high energy gamma rays (TeV range) were detected by us using an earlier, less sensitive, 
atmospheric Cerenkov telescope ( parraminana et al 1989 ; Brazier et al. 1990) and by another 
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group also using Cerenkov telescopes (Raubenheimer et al 1989). These observations, which were 
based on searches for pulsation, showed evidence for sporadic strong emission at a preferred orbital 



phase and for pulsations at, or in some cases near to, the pulsar period ( Lamb et al. 1977 ). We 
suggested that the emission may not be from a site close to the neutron star (Bowden et al. 1993). 
A recent paper gives evidence, from a re-analysis of archival data, for gamma ray emission from a 



site in the the accretion disc trailing the neutron star by 70° ( Raubenheimer fc Smit 1997 ). 



A GeV gamma ray outburst has been detected by the EGRET telescope on board CGRO 
( Vcstrand, Sreekumar &: Mori 1997a| ) at 100 MeV with evidence for pulsation at the pulsar period 
with an average luminosity ~ 5 x 10 36 erg s , a significant fraction of the total accretion luminosity 
which can vary from ~ 4 x 10 37 erg s _1 ( Finger, Wilson fc Fishman 1993] ) to 
(|White, Swank & Holt 1983| ). 



1.2 x 10 38 erg s" 1 



We report here our first observations with an imaging telescope of much higher sensistivity 
and a lower energy threshold than previously employed on this source. 



2. The Mark 6 Telescope 

The Durham Mark 6 atmospheric Cerenkov telescope has been described in detail elsewhere 
([Armstrong et al. 1997 ). It comprises three parabolic flux collectors of 7 m diameter and aperture 



f/1.0, mounted on a single alt-azimuth platform. At the focus of the central dish is a thermostat- 
ically controlled camera of 3° diameter comprising 109 photomultipliers in a close-packed array 
with hexagonal reflective cones to reduce light lost between photocathodes. The other two identi- 
cal parabolic collectors each have a temperature-controlled triggering camera containing 19 larger 
photomultipliers which cover an area similar to the central imaging camera. Events are recorded 
when there are responses in all three cameras which correlate in position at the foci (within 0.5° 
and in time to within 10 ns). This type of triggering gives this telescope a unique ability to detect 
low energy gamma rays from the ground with immunity from triggering by local muons, which are 
a significant limiting factor for low energy detection in single-mirror telescopes. The presence of 
independent noise in each of the triggered channels of the Mark 6 telescope results in a threshold 
which is not sharply defined. The probability of gamma ray detection reduces slowly from ~ 40% 
at an energy of 250 GeV to ~ 1% at an energy of 100 GeV. 

Comprehensive real-time calibration of the photomultipliers and digitizing electronics is per- 
formed by triggering the telescope at random times using (a) a nitrogen laser /plastic scintilla- 
tor/optical fibre guide/opal diffuser system to simulate Cerenkov flashes and enable flat-fielding, 
and (b) pulses injected after the coincidence stage to provide random samples of the background 
noise. The latter enables a correct pedestal or zero to be assigned to each fast ADC to an accuracy 
of ±0.1 digital counts. The laser trigger enables the gain of each photomultiplier to be measured 
relative to the others, to an accuracy of about 1% for each 15-minute data segment. Each type of 
random trigger is generated at an average rate of 50 per minute. 
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Simultaneous and continuous monitoring of the atmospheric clarity is carried out using an 
infra-red radiometer monitoring the source region and an axial optical CCD camera which enables 
the position and magnitude of guide stars to be measured. The digital data record for each event 
contains: (1) the digitized charge from each photomultiplier, (2) the time of each event, to an 
accuracy of 1 fis, using a GPS-moderated rubidium oscillator, and (3) the telescope pointing position 
to 1 arc minute. Every minute a record is made of the output of a comprehensive weather station, 
and of the temperatures inside the photomultiplier packages and the electronics, all important 
system voltages, and the anode currents and trigger rates of all photomultipliers. 



3. Observations 

Centaurus X-3 was observed using the Durham Mark 6 Cerenkov telescope for typically 3 hours 
per night, from 0900 UTC on 1997 March 1, 3 and 4 (corresponding to orbital phases (relative to 
X-ray eclipse) of 0.78, 0.74 and 0.22, respectively) and from 1030 UTC on 1997 June 1, 2, 4, 5 and 7 
(corresponding to orbital phases 0.89, 0.37, 0.33, 0.81 and 0.77, respectively). The majority of our 
observations were targeted to be near orbital phase 0.8 with respect to X-ray eclipse at which TeV 
emission had been detected earlier ( [Brazier et al. 1989 ). Data were taken in 15-minute segments 



alternately on- and off-source. The total cosmic ray background counting rate of the telescope was 
~ 400 per minute. The on-source segments were taken with Centaurus X-3 near, but not exactly 
at the center of the camera. The source was usually 0.1° — 0.2° from the camera center. Off- 
source data were taken by observing a section of sky at the same declination as Centaurus X-3, but 
separated by 15 minutes in right ascension. All data segments started so that the off-source and on- 
source segments possessed identical azimuth/zenith profiles. At the declination of Centaurus X-3 
(—60.6°), the center of the camera was pointed during the off-source segments at a location which 
was 1.84° from the source position. 

Data segments were accepted for analysis only if the sky was clear and stable and the total 
counting rates in each on/off-source pair were consistent at the 2.5 a level. After these requirements 
were satisfied a data set containing a total of 20 hours, equally divided between on-source and off- 
source observations, was analysed further. 



4. Analysis 

Routine reduction and analysis of accepted data comprises: 

1. calibration of gains and pedestals of all 147 photomultipliers and ADC electronics within a 
15-minute segment, using the embedded laser and false coincidence events, 

2. software padding of the data ( [Fegan 1997j ) to equalize the on-source and off-source photomul- 
tiplier noise, 
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3. identification of the precise location of the source in the camera's field of view for each event, 
using the axial CCD camera, 

4. a calculation of the spatial moments of the shower image relative to the source position, for 
each event, 

5. a rejection of events containing images which would be unlikely to be produced by gamma 
rays. 



The image parameters used for selection of gamma rays are found using the established 'image 
and border' technique (Fegan 1997). These parameters are: BRIGHTNESS, the total number of 



digital counts in the image, DISTANCE of the image centroid from the source position, ECCEN- 
TRICITY of the image (width/length), ALPHA, the angle by which the image's long axis misses 
the source and WIDTH, the RMS spread of the image along the minor axis. 

In addition to these parameters, we are able to exploit another property of gamma ray shower 
images: the relative freedom from fluctuations in the images from gamma ray showers in separated 
detectors. We had demonstrated this in the enhancement of the gamma ray signal from AE Aquarii 
using two well-separated (100 m) telescopes flChadwick et al. 1996j ). In the Mark 6 detector the left 
and right triggering cameras are 15 m apart and can accurately measure similar parameters. These 
parameters depend on only the lower moments of the images: brightness and centroid position in 



the focal plane. Of those developed in Chadwick et al. 1996, the parameter -D m i ss is least affected 



by the relatively small separation of the detectors; we denote the single-telescope analog of this 
parameter by -Ddist- This measures the angular separation in the focal planes of the centroids of 
the two images. The sensitivity of this measure to cascade development has been verified by Monte 
Carlo simulations of showers and by observing the variation in the distribution of the parameter 
with zenith angle, using off-source data. The expected variation of the parallax due to shower 
development changes with zenith angle is clearly seen. Observations of gamma rays from PSR 
B1706-44 confirm the value of the Q-factorL] for this technique to be ~ 1.4, in agreement with 
simulation results ( Chadwick et al. 1997a| ; |Chadwick et al. 1997b| ). 



5. Results 



Before background rejection, events are removed from the data which have images of very low 
brightness or whose centroids fall very close to the camera edge, more than 1.1° from the center of 
the camera. After these have been removed there are 64713 events in the on-source data set and 
65171 in the off, a source deficit of 1.3 a. The parameter limits were set as follows: 



defined as 



fraction of gamma rays retained 
^fraction of background retained 
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Fig. 1. — The difference of the distributions in ALPHA for on-source data and off-source data. 

1. 800 > BRIGHTNESS > 20000 total digital counts (corresponding to £ 400 GeV), 

2. 0.3 < ECCENTRICITY < 0.8, required to enable ALPHA to be estimated without large 
errors, 

3. 0.35° < DISTANCE < 0.75°, required to allow gamma ray shower images to be elongated 
(lower limit) and to reduce effects due to the edge of the camera (upper limit), 

4. D dist < 0.07°, 

5. WIDTH < 0.28°. 



After these selections, there were 4429 events on-source and 4087 off-source, an on-source excess 
of significance 3.7 a. There has been no normalisation of any of the on- and off-source parameter 
distributions. The distribution of ALPHA shown in Figure |l] gives clear evidence of the excess 
being at small values of ALPHA were the excess is expected from a true gamma ray source (see 
e.g. Fegan 1997 ). For ALPHA less than 30°, the on-source and off-source data sets contain 1546 
and 1208 events respectively, a 6.4 a excess: 338 ± 52 events. 

A false-source analysis has been performed by re-analysing the on- and off-source data using a 
matrix of trial source positions, in celestial coordinates. At each trial source position in the matrix, 
DISTANCE and ALPHA are recalculated for each event and a final selection of events is made 
with the DISTANCE selection as above and ALPHA < 30°. The excess events for each trial source 
position are shown in Figure The central point in the matrix is the position of Centaurus X-3 
and is not the center of the camera, which is typically off-source by 0.1° — 0.2°. The width of the 
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Fig. 2. — False source analysis: excess events with ALPHA less than 30°. The contours are spaced 
at 0.6 a intervals, with the grey scale being such that black indicates a probability of > 6a for a 
gamma ray originating from that direction. 

central peak is determined by the maximum value of ALPHA in the selection and the maximum 
and minimum values selected for DISTANCE. A very strong source would allow tighter selections 
and a narrower peak. 



5.1. Flux 

The 338 excess gamma rays were detected in 10 hours of on-source observation. The selection 
procedure retained ~ 50% of the original gamma ray events. The collecting area has been estimated 



using Monte-Carlo simulations ( Armstrong et al. 1997 ) and is ~ 10 5 m 2 at the energies and zenith 



angles at which Cen X-3 was observed. Using this estimate and assuming a differential flux index 
of —1.8, the gamma ray flux above 400 GeV was (2 ± 0.3) x 10 -11 cm -2 s . The flux would be 
~ 10% higher for an assumed spectral index of —2.6. The main uncertainty is not in the collecting 
area, which is dominated by geometric effects and the height of origin of the electromagnetic 
shower, but by the variation of triggering probability with energy. This may be estimated from 
simulations and checked by observations, but still may be in error by ~ 50%. Our result is shown 
in Figure || where it is compared with previous results, including the EGRET spectrum obtained 
from a recently reported burst of GeV gamma rays flVestrand, Sreekumar fc Mori 1997a| ). The 
burst occured during the week starting 1994 October 18 and had an integral flux above 0.1 GeV of 
9.2 ± 2.3 x 10~ 7 cm" 2 s _1 with a best fit differential index of -1.81 ± 0.37 in the range 0.1 to 10 



-7- 



E 
o 



X 

3 



CD 
O 




-2 -1 1 2 3 4 5 

LOG (gamma ray energy) (GeV) 

Fig. 3. — Gamma ray integral number flux spectrum. The solid line is the EGRET 1994 October 
burst ( |Vestrand, Sreekumar fe Mori 1997a| ). Key: • - this work, □ - upper limit from January 
1994 EGRET data ( |Vestrand, Sreekumar & Mori 1997a| ), ■ - |North et al. 1990| , o - |Yoshi et al. 



1996 



GeV. Figure || also shows a 2 a upper limit form an earlier EGRET exposure in 1994 January. 

Our data have been examined for short-term variations. The detected rate of gamma rays 
per hour is shown for the observations on eight nights in Figure ||. The night-by-night detected 
gamma ray rates show no significant deviation from constancy (x 2 = 9 for 8 degrees of freedom). 
We have also tested to see if there is any significant difference between the detected count rates for 
observations clustered around phase 0.25 and phase 0.75. We find that there is a 1.7<r difference to 
which we ascribe no significance. 



5.2. Periodicity 

A data set was made from on-source events selected to have ALPHA < 30° and the exact 
event times were reduced to the Solar System Barycenter using the DE200 ephemeris and then to 
the position of the neutron star using the current orbital ephemeris employed by BATSE ( [Bildsteri 
et al. 1997|) . 

Data from individual nights have been searched for periodicity using the Rayleigh test. The 
period search range was over one independent Fourier interval about the contemporary X-ray 
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Fig. 4. — Rate of detected gamma rays per hour for each night's observation. Data taken at orbital 
phases with respect to mid X-ray eclipse: o ~0.25, • ~0.75. 



period taken from the BATSE results^. The most significant Rayleigh power, at the X-ray period 
and corresponding to a chance probability of 3 x 10 -2 after correction for period searching, was for 
data recorded on 1997 March 3. One such occurrence is expectation for the number of observations 
made. The 3 a upper limit to the gamma pulsed fraction at the fundamental may be derived from 
the corresponding Rayleigh probability (= exp (-NR 2 )), given R, the length of the Rayleigh vector 



(Mardia 1972) and N, the number of events. 



The spin period history of Centaurus X-3 ( Bildsten et al. 1997 ) does not allow a search 
for coherent pulsations involving data sets of long duration. The data set has been split into 
two for coherent analysis, comprising the observations in 1997 March and 1997 June, and each 
was searched separately over a period range correponding to one independent Fourier interval. 
The X-ray period of Centaurus X-3 during each gamma ray observation has been taken from the 
contemporary BATSE results. There is no significant Rayleigh power at or close to the expected 
period or half-period in either of the the two datasets. The data sets from March and June had 
N = 806 and N = 740, leading to a 3<7 upper limit for the fraction of gamma rays which were 
periodic of 9.0% and 9.8% respectively, for either the X-ray period or the half-period. 

Following the earlier suggestions that gamma rays may originate elsewhere in the Centaurus 
X-3 system (Bowden et al. 1993, Raubenheimer & Smit 1997) we have analysed our data to search 
for evidence of pulsed emission at a site other than the neutron star. In this analysis, a matrix of 
points in the co-rotating frame of the system was set up and the times of the selected events were 
corrected for the orbital motion of each point. This analysis did not provide evidence of pulsed 



2 Available on the web at http:/ /www.batse. msfc.nasa.gov/data/pulsar. The parameters for Cen X-3 are based on 
observations with BATSE taken by Wilson et al. 
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emission at either the period or half-period from any region. 



6. Conclusion 

Centaurus X-3 is an extremely variable source of X-rays, varying on time scales from days to 
weeks. The periodic outburst reported by EGRET and the earlier failure to detect gamma rays 
show that this source varies in intensity by at least a factor of five at GeV energies over an interval 
of ten months ( Vestrand, Sreekumar fc Mori 1997a| ). The EGRET GeV outburst, wholly pulsed at 



half the X-ray period (Vestrand, Sreekumar &: Mori 1997b), was observed at a time of slow-down 



in the X-ray period when the neutron star was not accreting mass rapidly. Our 400 GeV unpulsed 
flux was observed at a time when the BATSE data indicated a period of spin-up (March data) and 
a transition between an interval of spin-down and an interval of spin-up (June data). 

Our results are the first measurements of the time-averaged very high energy (400 GeV) 
gamma ray flux from Centaurus X-3. The observation method adopted by the earlier, non-imaging 
Gerenkov telescopes could not detect such a continuous flux — a strong outburst of a pulsed signal 
was necessary to enable a detection. These earlier detections were confined to one night in a sample 



of 95 nights' data (Bowden et al. 1993) and one night from 59 nights' data (Raubenheimer &: Smit 



1997). The current detection of a lower level of continuous emission, stable over eight nights of 



observation, is not in conflict with these earlier results. 

The measurements summarised in Figure |3| suggest that Centaurus X-3 has a continuous 
unpulsed emission and sporadic outbursts of pulsed emission over a wide range of gamma ray 
energies. Models for gamma ray emission predict significant absorption from ~ 1 GeV to ~ 100 
GeV due to gamma-thermal X-ray interactions in the disc ( Bednarek 1993] ) . Future reductions in 



energy threshold of Cerenkov telescopes may enable the spectral variation in this important region 
to be measured. 

The earlier measurements of pulsed VHE gamma ray emission suggested that the pulsed emis- 



sion came from a localised area trailing the neutron star (Bowden et al. 1995; Raubenheimer fc 



Smit 1997| ). In contrast, the present detection of unpulsed emission indicates that this radiation 



must originate in an extended volume of the system. 
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